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We present a comparative assessment of the features of inelastic atom-surface scattering spectra 
that are produced by several different forms of linear and nonlinear phonon coupling to the projectile 
atom. Starting from a simple theoretical model of atom-surface scattering and employing several 
recently developed exact numerical and approximate analytical methods we calculate and compare 
the scattering probabilities ensuing from each form of interaction and from each calculational scheme. 
This enables us to demonstrate that in the regime of thermal energy atom scattering from surfaces 
the dominant contributions to the zero-, one- and multi-phonon excitation probabilities obeying 
unitarity arise from linear coupling treated to all orders in the interaction. 
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Theoretical interpretations of inelastic atom-surface scattering experiments carried out in the past two decades to 
investigate vibrational properties of the various surfaces [1] have all been based on several assumptions concerning 
atom-surface interactions that have never been critically examined or discussed. The aim of the present Letter is to 
remedy this situation in one of its so far neglected aspects, viz. the lack of a detailed assessment of the effects that 
arise from different forms of atom coupling to quantized surface vibrations or phonons. We perform this task using 
a simple paradigmatic, yet realistic enough quantum model comprising all the essentials of single- and multi-phonon 
atom-surface scattering. 

The first assumption adopted by the majority of the authors in the development of inelastic atom-surface scattering 
theories reviewed in Refs. [2-8] has been the linear coupling of scattered atoms to vibrational displacements of surface 
atoms (i.e. linear coupling to the phonon field), and only few authors have resorted to descriptions based on nonlinear 
coupling [9,10]. The absence of a quantification of the effects originating from the different forms of coupling is rather 
surprising because they should manifest themselves in all the components of the scattering spectra, viz. in (a) elastic 
scattering probability given by the so-called Debye- Waller factor (DWF), (b) one-phonon scattering probability, and 
(c) multi-phonon scattering probabilities. 

The second assumption concerns the treatment of one-phonon excitation probabilities. Here, a rather general 
consensus prevailed that these should be calculated in the limit of linear atom-phonon coupling and in the Born 
approximation using distorted waves accounting for reflections of the atom from the surface to all orders in the 
static atom-surface potential [5,8]. However, as the distorted wave Born approximation (DWBA) yields non-unitary 
scattering amplitudes it can provide reliable estimates only for the relative magnitudes of single quantum excitation 
probabilities in the regime in which the multi-quantum excitations are negligible. In order to obtain accurate one- 
phonon excitation probabilities for comparison with the zero- and multi-quantum ones (which are all sensitive to the 
form of atom-phonon coupling) and with the experimental data, the DWBA results should be unitarised to comply 
with the optical theorem. 

The most striking differences arising from the assumptions of cither linear or nonlinear atom-phonon coupling are 
expected in the calculated multiphonon scattering probabilities. Recent He atom scattering (HAS) experiments have 
successfully revealed and disentangled the one- and multi-quantum features in the atom-surface scattering spectra 
for phonons with Einstein [11,12] and acoustic [13,14] type of dispersion. Hence, the resolution of the origin of 
multiquantum features in the measured spectra as being due cither to (i) multiple single phonon scattering events 
describable in terms of linear coupling, or (ii) multiquantum excitations that arise from nonlinear coupling in each 
order of perturbation expansion, emerges as a task of utmost interpretative importance. 

One of the important conditions for a conclusive testing of the roles that linear and nonlinear coupling may play 
in the various theoretical descriptions of atom scattering from surface phonons is that they all ought to be applicable 
to a model which realistically describes the essentials of a particular scattering system. Another condition is that a 
clear disentanglement of the single- and multi-quantum features is possible for comparison of the results of different 
theories with the available experimental evidence. Experimental systems that enable the easiest identification of and 
distinction between the various one- and multi-phonon excitations in atom-surface collisions are those in which atoms 
are inelastically scattered by dispersionless or Einstein phonons. In this case the inelastic scattering spectra display 
a series of well separated equidistant peaks at the multiples of Einstein phonon energy ±fi,wo away from the elastic 
line. The weight of the n-th peak measures the probability of excitation of n phonons in the course of scattering. 
A prototype collision system of this kind is provided by HAS from monolayers of Xe atoms adsorbed on Cu(lll) 
and Cu(001) surfaces which support vertically polarized Einstein modes [11]. Other systems with similar properties 
have also been reported in the literature [10,12]. The model Hamiltonians required for a full quantum mechanical 
description of these systems have been presented elsewhere [2,11]. However, for the present purpose all the relevant 
features of collision dynamics are included in the effective one-dimensional Hamiltonian H that describes motion of 
a projectile nonlincarly coupled to a surface oscillator of frequency ui - Denoting by z (Z), p (P) and m (M) the 
coordinate, momentum and mass of the scattered particle (surface oscillator), respectively, we can write 



Here, the elastic particle motion in the model static surface potential of the Born-Mayer form, U(z) = De~ az , is 
described by distorted waves [15] that are eigenstates of the particle Hamiltonian H*(p,z) = p 2 /2m + U{z). By 
introducing a purely repulsive U (z) we avoid the occurrence of resonant scattering effects associated with the bound 



p^_ 

~ 2m + 2M 





+ De- az (e aZ - 1) 



(1) 



2 



states of the surface potential well [16-18] that are unimportant for the present study. The inelastic scattering is 
governed by the interaction potential 

V(z,Z) = De- az (e aZ - 1) 
Z 2 a 2 

= ZaDe- az + — ^— De- az + ... 

= ZV 1 (z) + Z 2 V 2 (z) + ..., (2) 

that contains powers of oscillator displacements Z to all orders. The main merit of the thus defined model lies 
in its amenability to exact numerical [9,19] and approximate analytical treatments [2] in the calculations of single- 
and multi-phonon scattering probabilities. This enables us to consistently investigate as how the different forms of 
atom-phonon coupling affect the scattering spectra. 

In this Letter the exact numerical treatments of elastic and inelastic atom scattering probabilities are carried out 
by employing the coupled channels method [19] (CC) to solve the Schrodingcr equation of particle motion perturbed 
by three types of interactions embodied in expression (2): the linear coupling term ZVi(z), the sum of linear and 
quadratic coupling terms ZVi(z) + Z 2 V 2 {z), and the full potential V(z,Z). Tractable analytical solutions for the 
same model are obtained for the cases of linear coupling ZV\(z) treated to all orders within the distorted wave 
exponentiated Born approximation (EBA) formalism [2] based on cumulant expansion, and for the full coupling 
V(z, Z) treated in the DWBA as regards the particle motion. With the exception of the latter approximation, all 
other mentioned treatments yield unitary results in that the sum of all calculated scattering probabilities equals unity. 
Separate treatment of the sum of linear and quadratic coupling terms is particularly instructive because the quadratic 
term gives rise to both inelastic and elastic scattering processes, viz. the simultaneous two phonon excitations and 
clastic scattering by the static potential V^z) renormalised by emission and reabsorption of a virtual phonon, an effect 
that is not present in linear coupling theories. These features are diagrammatically illustrated in Fig. 1 in which the 
virtual phonon excitation is represented by a closed phonon loop. Diagrammatic perturbation expansion in V(z, Z) 
shows that nonlinear coupling to all orders gives rise to the appearance of arbitrary number of closed phonon loops 
in any single- or multi-phonon vertex [20]. This leads to a Holstein type of renormalisation [21] of inelastic particle 
scattering matrix elements by the factor exp[— (crao) 2 /2] where Uq is the mean square displacement of the oscillator. 

The calculations have been carried out for the elastic scattering probability Po, the one-phonon scattering probability 
Pi, and the two-phonon scattering probability Pi for the types of couplings and within the theoretical approaches 
described above by using the parameters representative of He and Ne atom scattering from adsorbed Xe atoms. Figure 
2 compares the behaviour of P n (n — 0, 1, 2) calculated for linear, linear and quadratic, and full nonlinear couplings as 
the functions of scattered particle incoming energy E i} for He atom scattering from Xe atoms adsorbed on a cold Cu 
surface. The generalisation to finite surface temperatures is straightforward [9,16] but does not bring in new effects 
concerning the various types of atom-phonon coupling. 

The important message that can be deduced from the behaviour of elastic scattering probability P or the Debye- 
Waller factor depicted in Fig. 2 is an excellent agreement between the results of exact CC calculations for the sum of 
linear and quadratic and full nonlinear coupling in the whole range of the projectile incoming energy. This means that 
the nonlinear coupling effects are described to a high degree of accuracy already by the sum of linear and quadratic 
coupling terms. The EBA formalism [2] that is based on approximate treatment of linear coupling to all orders is also 
in excellent agreement with the exact numerical results. 

Almost identical conclusions pertain to the behaviour of the one-phonon scattering probabilities Pi in Fig. 2. As 
expected, the exact linear coupling result slightly deviates from the exact nonlinear coupling results as it does not 
pick up all the processes that are present when the nonlinear coupling terms are switched on. Here we also show the 
one-phonon scattering probability calculated in the DWBA applied to the full nonlinear coupling interaction V(z, Z). 
For very low incoming energies near the one-phonon excitation threshold E^ = Kloq = 2.7 meV, i.e. in the one-phonon 
scattering limit in which linear coupling yields the dominant contribution, the scattering probability calculated in the 
DWBA with full nonlinear coupling is in a good agreement with the other ones. However, with the increase of P, 
the deviations from the exact result soon become large and signify the break down of this approximation due to the 
nonunitary treatment of all the scattering events induced by V{z, Z). 

Very similar trends are observed in the behaviour of two-phonon scattering probabilities (P2 in Fig. 2). Again, the 
EBA results are in a very good agreement with the exact ones and hence this formalism proves to be a very reliable 
approximate method for treating multiphonon scattering in the quantum regime. This property allows us to use the 
analytical form of the EBA scattering probabilities [2] to demonstrate that the maxima of Pi(Pj) and p2(Pi) appear 
approximately at the values of Ei for which Pi = P and P 2 = Pi, respectively, and that Pi (Pi) is largely given 
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by P\{Ei) — n(Ei) exp[— n{Ei)\ where n(Ei) — — In Po(Ei) is the mean number of phonons excited in the scattering 
event [2]. 

Qualitatively very similar results are retrieved for Nc atom scattering treated in the same model (see Fig. 3). 
Likewise the case of He scattering, all the unitary treatments yield quantitatively similar results which follow the 
general trends present in Fig. 2. The only strong discrepancy between exact and approximate results appears for 
the scattering probability calculated in the DWBA with nonlinear coupling, that again completely fails to reproduce 
exact results. The centre and lower panels also demonstrate as how due to the larger projectile mass the maxima of 
all the unitary one-phonon and two- phonon scattering probabilities shift to lower Ei than in the case of lighter He 
atoms. This is readily explained by the scaling properties of n(Ei) which in the quasi-elastic limit of the EBA, here 
valid for Ei > 10 meV, is very accurately reproduced by the expression n(Ei) = 8uQmEi/h 2 that is independent 
of both a and D. This causes exponentially faster attenuation of the Debye- Waller factor and all phonon excitation 
probabilities for larger projectile masses and fixed incoming energies. 

The most striking feature of the results shown in Figs. 2 and 3 is a complete failure of the DWBA with full 
nonlinear coupling to describe the two-phonon scattering probabilities. In view of this and the earlier established 
excellent agreement between the measured and calculated EBA scattering probabilities for several prototype systems 
[2,11], the most important conclusion that can be drawn from Figs. 2 and 3 is that in inelastic low energy atom-surface 
scattering the main contribution to the two- and higher order phonon scattering probabilities comes from successive 
one-phonon scattering events and not from the processes of simultaneous multiple phonon excitations from a single 
interaction vertex that arise in nonlinear coupling models. This clearly demonstrates that unlike the case of inelastic 
neutron scattering from crystal lattices, the application of DWBA combined with nonlinear coupling to interpret the 
multi-phonon features in thermal energy atom scattering spectra leads to erroneous results and therefore should be 
avoided. 

In the present case of a one-dimensional scattering model the CC-computations of the scattering probabilities are 
feasible for both the linear and nonlinear forms of coupling. Therefore, in such cases it is advantageous to immediately 
resort to the full nonlinear coupling in the CC-calculations. This ceases to be the case with the treatments of nonlinear 
coupling by other methods that yield scattering probabilities fulfilling the unitarity requirement. Hence, the finding 
that the EBA results arc in a very good agreement with the exact ones becomes particularly important for the 
treatment of full three-dimensional scattering problems involving dispersive phonons for which the CC-calculations 
may not be feasible yet. In view of this, the present results affirm the EBA as a powerful tool for reliable interpretation 
of low energy atom scattering from surface phonons [2] . 

In summary, we have investigated the effects of different forms of linear and nonlinear atom-phonon couplings on 
inelastic atom-surface scattering. By comparing the results of several exact numerical and approximate analytical 
treatments of a simple model of atom scattering from surface vibrations we have demonstrated that in the regime of 
thermal energy atom scattering the linear coupling interaction treated to all orders gives a dominant contribution to 
the scattering probabilities. 
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FIG. 1. Diagrams for the diagonal components of the scattering matrix whose imaginary parts yield the total probability 
for the processes induced by the depicted phonon emission and absorption vertices. Full and wavy lines denote the atom and 
one-phonon propagators, respectively. Open circles, full circles and full squares denote vertices proportional to (uoa), (uoa) 2 , 
and (moq) 4 , respectively, (a) One-phonon scattering process induced by the linear coupling term ZVi(z) treated to second 
order, (b) A two-phonon scattering process encompassing one "static vertex" Vi(z) renormalized by a phonon closed loop. 
This diagram is induced by the quadratic term Z 2 V2(z) treated to third order, (c) A two-phonon scattering process with 
renormalized two-phonon vertices that is induced by the quartic coupling Z 4 V4(z) treated to second order. 

FIG. 2. Scattering probabilities for elastic (upper panel), one-phonon inelastic (centre panel) and two-phonon inelastic (lower 
panel) He atom scattering from a vibrating Xe atom on a cold surface, calculated for the three types of interactions within the 
scattering formalisms discussed in the text and shown as functions of the incoming He atom energy Ei. Solid line: CC result for 
full V(z,Z); dotted line: CC result for ZVi(z); dashed line: CC result for ZVi{z) + Z 2 V 2 (z); dashed-dotted line: EBA result; 
chained squares: DWBA result for full V(z,Z). Note that P = 1 and Pi = below the one-phonon excitation threshold at 
huo, and P2 — below the two-phonon excitation threshold at 2fioJo- 

FIG. 3. Same as in Fig. 2 but for Ne atom scattering from a Xe atom on a cold Cu surface. Comparison with Fig. 2 shows 
how the scattering probabilities scale with the projectile mass when all other parameters are fixed. 
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